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ABSTRACT 

The modified version of a bullet model for gamma ray bursts is studied. The central 
engine of the source produces multiple sub-jets that are contained within a cone. 
The emission of photons in the source frame of a sub-jet either takes part in an 
infinitesimally thin shell, or during its expansion for a finite time. The analysis of the 
observed profiles of GRBs taken by BATSE leads us to the conclusion that the latter 
possibility is much more favored. We also study the statistical distribution of GRBs, 
in the context of their bimodality of durations, taking into account the detector's 
capability of observing the signal above a certain flux limit. The model with shells 
emitting for a finite time is able to reproduce only one class of bursts, short or long, 
depending on the adopted physical parameters. Therefore we suggest that the GRB 
bimodality is intrinsically connected with two separate classes of sources. 

Key words: gamma rays: bursts - gamma rays: theory 



1 INTRODUCTION 

The observed durations of Gamma-ray Bursts range from 
milliseconds to several hundreds of seconds, forming two 
distinct peaks. Thus two GRB classes are proposed: short 
(T 2 s) and long (T ^ 2 s) bursts (Kouvelietou et al. 
1993). (However, there have recently been claims of the 
third, intermediate duration peak; see Horvath et al. 2004). 
The coUapsar scenario (Woosley 1993; Paczyriski 1998) is 
commonly favored for the origin of a long GRB event, based 
on observations of afterglows and host galaxies of these 
events (see the review by Zhang & Meszaros 2004; Piran 
2004). In particular, in some long bursts the afterglow obser- 
vations clearly show the association with supernovae (Stanek 
et al. 2003) , thus confirming the coUapsar scenario. Further- 
more, the GRB positions inferred from the afterglow obser- 
vations are consistent with the GRBs being associated with 
the star forming regions in their host galaxies. 

Yet, the situation is still far from clear for the short 
events, and various possibilities have been discussed. 

An essentially different mechanism, such as a compact 
binary merger (Eichler et al. 1989; Paczyhski 1991; Narayan, 
Paczyhski & Piran 1992) was proposed. The duration of the 
central engine activity is in this case of the order of 2 sec- 
onds, appropriate to account for the fuefing time of a short 
GRB. For many years, no optical counterpart was observed 

* E-mail: agnes@camk.edu.pl 



for a short burst but the first detection of the X-ray after- 
glow by SWIFT satellite indicated a location 9.8 arcseconds 
from the center of the elliptical (type El) galaxy (Burrows 
et al. 2005) at a redshift of 0.2249 -I-/-0.0008 (Prohaska et 
al. 2005). In such a galaxy strong starburst activity is not 
expected, and no traces of the relative supernova were found 
in the optical data (e.g. Kosugi et al. 2005) which again may 
favor a merger scenario. The results and implications for the 
burst GRB 050509B were recently discussed in Gehrels et 
al. (2005). For the recent observations of the GRB 050709, 
the source can be associated with the star forming galaxy at 
z=0.16 (Fox et al. 2005; Villasenor et al. 2005), however with 
no supernova associacion in the optical lightcurve (Hjorth et 
al. 2005). Taking into account the distance, the luminosity 
of this GRB afterglow is about 3 orders of magnitude lower 
than for a typical long event, making the coalescing compact 
binary the most promising progenitor candidate (see also 
the discussion by Piro 2005). Furthermore, recent observa- 
tion of of GRB 050724 afterglow emission supports strongly 
the merger origin of this event (Berger et al. 2005). 

Also, the events of neutron star-neutron star (NS-NS) 
or neutron star-black hole (NS-BH) mergers should hap- 
pen quite often in the Universe (i.e. Bulik, Belczyhski & 
Kalogera 2003; Bulik, Gondek-Rosihska i Belczyhski 2004). 
These events should have given us a detectable sign of their 
presence, not only in the form of gravitational waves, but 
also in the gamma-rays. 

On the other hand, apparently the main argument in 
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favor of the merger scenario is that we cannot reject this 
possibihty. Therefore there have been attempts to unify one 
model for all the GRBs, and assign all the bursts to the 
most plausible scenario, i.e. coUapsar. One of the important 
arguments here is the observation that a short GRB may be 
similar to the first ~ 1 s of a long one (Ghirlanda, Ghisellini 
& Cellotti 2004). 

It has been recently suggested (Yamazaki et al. 2004) 
that the bimodal duration distribution of GRBs can be ex- 
plained in the frame of inhomogeneous model of a GRB 
jet consisting of multiple sub-jets, emitting for an infinites- 
imally short time. Since the sub-jets, or sub-shells, are not 
distributed uniformly within the main jet, but instead may 
have e.g. Gaussian (or power-law) distribution depending 
on the angular distance from the jet axis, the observer does 
not always detect the same number of sub-jets on his line 
of sight. In other words, the number of sub-jets, ris, seen by 
the observer depends on the viewing angle. 

The multiple sub-jet model proposed by Yamazaki et 
al. (2004) assumes that the duration of the observed burst 
does not essentially depend on the duration of activity of the 
central engine. Or, at least, the central engine of both short 
and long bursts may be of the same nature, i.e. a coUapsar, 
operating for, say, 30 seconds. It means that duration of ac- 
tivity of the central engine is only slightly reflected in the 
duration of the observed GRB event. Instead, the event du- 
ration depends mainly on the viewing angle. This is because 
the observer may either detect a large number of sub-jets 
(ws > 1), and in consequence the burst duration is long, or 
the number of sub-jets detected on the line of sight is tIs = 1 
and we observe a short burst. 

In this article we build a similar model of a non-uniform 
jet, that consists of a number of randomly distributed sub- 
jets. We study both the case of sub-jets emitting the gamma 
rays at an instantaneous time t = to, aa invoked by Ya- 
mazaki et al (2004), and the case of a more realistic assump- 
tion that the emission of gamma rays is not instantaneous 
in time, but lasts for a certain period during the expansion 
of the shell from i?min to i?max. We calculate the profiles of 
individual sub-bursts as well as the total lightcurves of the 
GRBs that are the sum of multiple sub-bursts. Next, we fit 
the profiles of the individual pulses using the phenomeno- 
logical prescription of Ryde & Svensson (2002) and check if 
these profiles correspond to these seen on the observational 
data from BATSE. The data were analyzed both by Ryde 
& Svensson (2002) in case of long duration events and by 
ourselves in case of short events. 

We also simulate the statistical distribution of durations 
of the GRBs. Firstly, we hold the assumption, that only the 
sub-jets present on the line of sight of the observer account 
for the GRB event. Secondly, we release this assumption, 
and include also the emission from the sub-jets detected off- 
axis. The statistic is calculated on the basis of the condition 
used for classifying the events as GRB, i.e. the flux limit 
(the detector's capability). 

In Section |5| we give the details of the model and as- 
sumptions. In Section l3.ll we present the resulting profiles 
of the bursts and sub- bursts and in Section [3.21 we compare 
these profiles with observations. In Section f3.3l we show the 
statistical distributions of the burst durations that result 
from the multiple sub-jet model. Finally, in Section |4] we 
discuss our results and give conclusions. 



2 MODEL 
2.1 Geometry 

We consider a coUimated jet in the form of a cone with the 
opening angle A^tot, which is present for the whole period 
of the duration of central engine activity, Tdur. The main 
jet consists then of a number A'sub of sub-jets, which have 
smaller opening angles A^sub and are launched at their char- 
acteristic departure times, extending until the time Tdur: 
< tdcp < Tdur- The position of each sub-jet is determined 
by two angles, (6i, tpi), such that the angular distribution of 
the direction 6i, i.e. the distance from the z axis of the main 
jet, is given by the Gaussian function: 

P(ei)dei = e-4fe'd6(i, for < A^tot - Ae,,b (1) 

where 6c < A6'tot is a parameter of this distribution. Inde- 
pendently, the angle ipi, measured in the x — y plane in the 
basis of the jet, is uniformly distributed between and 27r. 

We assume that all the subjets have the same Lorentz 
factors 7 and up to the radius ro, for a period of time to, 
expand in a photon-quiet phase. At the radius ro the photons 
are emitted instantaneously, i.e. the duration of the photon- 
active phase in the comoving frame is infinitesimally small. 
The sub-jet i is expanding for the time: 

to ~ tdcp + (2) 

where /3 — -^1 — 1/7^- Then, all the photons are emitted in 
an infinitesimally short time, at: 

= ^ (3) 

after the launch of a sub-jet. 

The observer's position is determined by his viewing 
angle (Sobs, ^ohs), as measured in the reference frame of the 
main jet. Again, we randomly choose these angles, without 
favoring any particular direction, and the angles are < 
Sobs < A^tot and < ^obs < 27r. The observed temporal 
structure of a flash of radiation produced by a sub-jet is 
determined by the time measured in the observer's frame, 
i.e. the arrival time of a photon: 

= (1 + ^) (to + 4.(1 -13 cos \)) (4) 

where we take into account that relativistic aberration of 
photons emitted at an angle A to the observer occurs only 
after the expansion time to. This is because the assumed 
radius of the emitting photosphere is large and cannot be 
neglected. Here the angle A is measured between the direc- 
tion of motion of the emitting plasma with respect to the 
observer. 

Because of the axial symmetry of the problem we choose 
an observer to be located at i^obs = in the reference frame 
of the jet. In order to determine this angle A and its depen- 
dence on the observer's position, location of the subjet, and 
direction of the photon, let us consider a spherical triangle 
located on the surface of the sphere of the radius ro. The 
vertexes of this triangle are determined by the directions 
of the jet axis, observer's axis, and sub-jet axis. Therefore 
the sides of the triangle are Sobs, Si, and Ai, while the angle 
opposite to Ai is i^i (see Figure 0. 

This implies the dependence: 

cos Ai = cos Si cos Sobs + sin Si sin Sobs cos(<^obs — ^i) (5) 
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Figure 1. Geometry of the problem. The jet has a half-opening 
angle A^tot- The observer is located at (ipobs = 0, Sobs) in the jet 
reference frame. Each subjet has a half-opening angle AS; and is 
located at {ipi,9i). The angle between the observer and the sub-jet 
axis, Ai, form a spherical triangle together with the angles Sobs 
and 9i. 

Now, the maximal and minimal angle A between the 
direction of motion of the emitting plasma with respect to 
the observer, is given by: 

Amax = + A^sub (6) 

Amin = mm(0; Ai - A9sub) 

Having calculated the angle Ai, we can determine the 
time Tobs at which the photon hits the detector. The du- 
ration of a single sub-burst is therefore given by the differ- 
ence between the ending and starting times in the observer's 
frame: 



an exact formula for the angle A but rather its upper and 
lower limits resulting from the opening angle of a sub-jet. 
The derivation of an exact expression would require a trans- 
formation between the two reference frames: jet's (in which 
the observer's position is given) and sub-jet's (in which the 
photons from a single sub-jet are produced). 

Therefore the simplified condition for the observer's po- 
sition to fall inside the sub-jet is: 

Ai < A^sub. (8) 
and also neglects the dependence on the angle ifi. 

2.2 Duration of a burst 

If the line of sight of the observer coincides with only one 
sub-jet, then the duration of this sub-jet is equivalent to the 
duration of the GRB. Alternatively, the observer may 'catch' 
a number of sub-jets. In this case the duration of the GRB 
is the sum of durations of all the sub-bursts that coincide 
with the line of sight: 

AT = Max{TU) ~ Min{T:,^,,) (9) 

The absolute maximal and minimal duration of a GRB 
event are estimated as follows. The duration of a single sub- 
burst for Ai > A6'sub is: 

cAt = ro(cos(Ai - A^^ub) - cos(Ai + A6'eub)) (10) 
= 2ro sin Ai sin ASgub 

while for Ai < A6'sub it is: 

cAt = ro(l-cos(Ai + A6'sub)) (11) 

(cf. Eq. 01. Therefore the shortest possible sub- bursts are 
these seen on-axis (Ai = 0), and: 

AWn^ (A^sub)' (12) 

c Z 

For our standard parameters, tq — 10^"* cm and A6'sub ~ 
0.02, equation (I12II gives Atmin = 0.66 s. 

The longest theoretically possible burst would be seen 
at the angle Ai = 7i"/2, and when the sub-bursts caught on 
the line of sight were so numerous that they would cover 
the whole jet opening angle: instead of A^sub we take AStot. 
In case of our parameters (AStot = 0.2, Tdop = 30s) this is 
about 1350 s. In practice, as discussed in the next section, 
the light aberration will effectively limit the observed surface 
to within a small angle around the line of sight. In such a 
case the duration of the longest burst will strongly depend 
on the central engine lifetime Tdep and on the sensitivity of 
the instrument used to observe the burst." 



riart = To + (1 + 2)^(1 - /3C0S A^in) (7) 
Tend = n + (1 + Z)g(l - /3C0S A„,ax) 

where T(j = (1 + z)4- 

In the above considerations, for simplicity we neglect 
the dependence of the angle A between photon's direction 
and observer's position on the angle ipohs- This is because we 
have random observers uniformly distributed between and 
27r, so for a sample with large enough statistics this will not 
influence the results. Also, please note that we did not give 



2.3 Sub-burst profiles 

The duration of the GRB event depends not only on the 
number of sub-jets that are observed, but also on the en- 
ergy flux that is measured by the detector. This is of par- 
ticular importance if we allow also for the sub-bursts that 
are observed off-axis, to contribute to the observed spec- 
trum. Because the energy flux observed from the relativistic 
jet sharply decreases for the viewing angles that are outside 
the jet cone, only a part, if at all, of such a sub-burst could 
be classified as a gamma ray burst. 
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The energy flux emitted from an expanding shell, as 
seen by the distant observer located on the jet axis, was 
studied in Fenimore at el. (1996). They considered both the 
case of an infinitesimally thin and extended shell. 

The photons emitted from a relativistically expanding 
shell, that arrive to the observer at the same time, originate 
from a prolate ellipsoid (Rees 1966). In order to calculate 
the flux in the detector rest frame, we have to integrate the 
emitted spectrum over the surface resulting from the cross- 
section of the ellipsoid with the sphere of the center in the 
detector and radius equal to the distance to the source: 

dF^{T) <x (j) f{u')S^dA (13) 

where S = [7(1 — ficosO)]^^ is the Doppler factor, fiv') is 
the emissivity in the comoving frame, and v' = {l + z)i'"f{l~ 
(3 cos 6). The surface differential is determined by: 

dA= ^A<i> sin ede (14) 

where 

A$ = TT if e< (A^eub - e^bs) or eobs = 0. 

/ cos A^sub — cos 61 cos 0' h. , , . 

A$ = 2arccos( ^) otherwise (15) 

smt^smt^obs 

The angle S^bs = ''^i is between the direction to the 
observer and the sub-jet axis. The angle 6 depends on time 
in the observer's frame, and here is how the observed pulse 
varies with time: 

e(r) = arccos(i--^-(T~ro)) (16) 
p I + z r 

We calculate the pulse proflle for Tstart < T < Tend, and 
in principle we should integrate from 6'min to Smax (where 
Omin = if ^obs < A6'sub or 6'niin = 6'obs " ^^sub Otherwise, 
and Smax = ^obs + ^^sub)- Howevcr, in case of an infinitesi- 
mally thin shell, we must keep r — ro and only one value of 
6 will satisfy the relation il6l In other words, the integral in 
equation 1131 becomes one-dimensional integral over the arc 
of the length Acp. 

2.4 Shell emission for a finite time in the 
comoving frame 

The model described above considered the emission that is 
infinitesimally short in the comoving frame, and the dura- 
tion of the observed pulse is a result of the purely geomet- 
rical transformation to the observer's frame. In reality, the 
emission from the expanding shell may not be instantaneous, 
but rather lasts for a certain time. In our model the du- 
ration of the emission is determined by the minimum and 
maximum radii of the expanding shell, Vmin and rmax- The 
starting and ending time of an individual sub-burst in the 
observer's frame is given by: 

TAart = (1 + ^)(tdcp + - /3C0S A„,i„)) (17) 

Tend = (1 + Z)(tdcp + -^(1 - 13 COS Amax)) 

where Amin and Amax are given by the equations |S| In or- 
der to obtain the time profile of each sub-burst, we have to 
calculate the integral in Eq. 1131 over 9 and A<j>, where the 
limits for 0{T) will satisfy the condition r-min < r < rmax- In 



this way we sum up in the observer's frame all the photons 
that are emitted in different moments in the expanding shell 
between rmin and Tmax but arrive to the detector at the same 
time T. 

2.5 Comparison with previous work 

The idea of the GRB jets consisting of multiple sub-jets was 
proposed by Heinz & Begelman (1999), and recently mod- 
eled in detail by Yamazaki et al. (2004) and Toma et al. 
(2005). In the latter work the authors successfully repro- 
duced the bimodal distribution of the GRB durations. Here 
we basically follow the same scheme of the multiple sub-jet 
model, however there are a few differences. 

(i) First, we calculate the angles between the emitting 
electrons and observer's direction, explicitly taking into ac- 
count the sub-jet azimuthal angle, ipi (see equation |S] and 
Fig.EJ. 

(ii) Secondly, we distinguish between the cases with and 
without the off-axis emission included in the total GRB 
pulse. This is particularly important, when determining the 
duration of a short burst and statistical fraction of short 
bursts in the simulation. 

(iii) Therefore, thirdly, we adopt a different way of se- 
lecting the pulses that are responsible for gamma ray signal. 
Instead of the spectral hardness, we adopt a criterion based 
on the sufficient signal with respect to the background fiux. 
The minimum flux is inversely proportional to the square 
root of Tgo, as implied by the BATSE characteristics (see 
Section 13.3.21 . 

(iv) Finally, we include in the model the possibility of 
a non-instantaneous emission, in the expanding shell that 
emits radiation from Rmin to -Rmax. 



3 RESULTS 

3.1 Pulse profiles 

We present here the shapes of the profiles for individual 
sub-bursts as well as the total GRB events that consist of 
multiple sub-jets. First, we show the results for the case of an 
infinitesimally thin emitting shell. In this case the emission 
in the source frame is instantaneous, and the duration of the 
event is determined only by the geometry. 

3.1.1 Profiles for the infinitesimally thin shell emission 

In Figure|5|we show the exemplary profiles of the pulses from 
the individual sub-jets, as seen by the observer. The sub-jet 
opening angle is A^aub = 0.02 and therefore the pulses that 
are seen inside this opening angle are much more energetic. 
The pulses seen off-axis are weaker, due to the strong de- 
pendence of the fiux on the Doppler factor and in turn on 
the viewing angle. This means that the maximum of the 
observed fiux is almost an order of magnitudes lower, as 
soon as the observer is located outside A^^ub (see e.g. cases 
of ^'obs = 0.018 and 6'^^^^ = 0.023). On the other hand, for 
larger viewing angles the pulse duration, as well as its start- 
ing time, substantially increase. For 6ohs = 0.007 a break in 
the light curve is caused by the conical geometry of the jet 
as discussed by Rhoads (1997). 
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Figure 2. The exemplary pulse profiles emitted from the indi- 
vidual sub-jets as seen by the observers located at various lines of 
sight. The angle 6»^^^ = 0.007,0.018,0.023 and 0.126 is indicated 
in each panel. The sub-jet opening angle is always AS^ub = 0.02 
and the Lorentz factor is 7 = 100. The flux is given with arbitrary 
normalization and essentially does not depend on frequency (flat 
spectrum) . 



Figure 3. Profiles of the GRBs for various observing angles with 
respect to the main jet axis: eobs=0.008, 0.054, 0.082 and 0.123. 
The lightcurves are plotted in 4 energies: 1.41 (solid line), 112 
(short dashed line), 482 (long-dashed line) and 2070 keV (dotted 
line). The bursts contain the total emission (on and off-axis) from 
the sub-jets, emitting at ro = 10^^ cm, that satisfy the condition 
Eq.|^(see text). 



In Figure 121 we show the total profiles of GRB, which 
are the result of multiple pulses seen by the observer. The 
plots axe labeled here with the observer angle with respect to 
the axis of the jet. While looking at small Sobs, the observer 
will detect a large number of sub-pulses that overlap, and 
give rise to a long, variable GRB. In contrast, the observer 
looking at larger angle 6lobs, can detect only a few sub-pulses. 
These form either a set of separate narrow spikes, or are 
connected with each other by the emission from these sub- 
bursts that were seen only off-axis (see Section 13.3.21 and 
condition Eg. 1211 for the detection of the off-axis emission). 

The spectrum in the comoving frame is defined using 
the broken power law function (see also the formula in Band 
et al. 1993): 

= i'^'/'^'o)^'^"' for v' < ly'o 
{i^'M'^" for i^'>v'o (18) 

In the comoving frame the break frequency is assumed v'q — 
2keV, and the spectral indices are a = —1 and (3 = —3 (the 
averaged values of the observed spectral indices were given 
e.g. by Pendleton et al., 1994; Preece et al. 1998, 2000). The 
lightcurves in Fig.j^are plotted in 4 energies: 1.41,112, 482 
and 2070 keV. 

3.1.2 Emission from the thick expanding shell 

Now, we show the results for the case of a thick expanding 
shell. In this case the duration of the individual pulse is de- 
termined not only by the geometry but also by the duration 
of the emission in the source frame. 



In Figure 21 we plot several examples of the individual 
sub-bursts, for various observing angles with respect to the 
sub-jet axis: S^bs = 0.004, 0.017, 0.024 and 0.127. 

In Figure |^ we plot the total GRBs, consisting of 
multiple sub-pulses, for several exemplary values of the 
observers angle (with respect to the main jet axis): 
e^bs = 0.009,0.044,0.096 and 0.128. The Figure shows the 
lightcurves in 4 energies: 1.41, 112, 482 and 2070 keV. The 
jet opening angle is A6'tot = 0.2. For the bursts that are 
observed close to the edge of the cone (large Sobs), the prob- 
ability of detecting a sub-jet on-axis is very small. These 
bursts consist mostly of the off-axis emission from the sub- 
jets, typically either with one-two strong, on-axis pulses 
(case of 6'obs = 0.096) or with none on-axis pulse (case of 
Sobs = 0.128). These bursts have the longest total durations- 
and last for several hundreds of seconds. The bursts that are 
observed closer to the main jet axis, are stronger and more 
variable, and their durations typically range up to ^ 100 
seconds. 



3.2 Fitting the pulse profiles 

The observed profiles of GRBs can be fitted to the exponen- 
tial function in the rise phase: 

N{t) (X exp(mt) (19) 

and to the power law shape in the decay phase (Ryde & 
Svensson 2002): 

iV(t) oc (l-ft/r)" (20) 
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Figure 4. Profiles of the individual sub-bursts for various ob- 
serving angles with respect to the sub-jet axis: 6'^j^_,=0.004, 0.017, 
0.024 and 0.127. Thick expanding shell approach. 



Figure 5. The same as in Fig. |3 the observing angles are: 
eobs=0.009, 0.044, 0.096 and 0.128. The bursts contain the total 
emission (on and off-axis) from the sub-jets, from i?min = lO'^* 
cm to -Rmax = 2 X lO'^* cm. 



where m,n and r are the free parameters. As shown by 
Ryde & Svensson (2002) for a sample of BATSE data, the 
distribution of the index n, characteristic for the pulse decay 
phase, peaks twice: around n = 1 and n = 3. However, 
one should keep in mind that this observational result has 
been obtained using a limited statistics of ~ 20 long GRBs 
detected by BATSE. Here we fit the pulse profiles that are 
obtained from the sub-jet model. 

First, we fit the pulses resulting from the assumed delta- 
function emission {t = to, r — ro) in the comoving frame. 
In the Figure |H| we show the distribution of indices for the 
sub-pulses that fitted to this pulse profile, and for which 
the maximum gamma-ray flux was substantially large (this 
corresponds roughly to the observing angle O'^y^^ > 0.03). 
The indices in our distribution are in the range n = 4 — 6, 
while neither n = 1 nor n = 3 showed up, contrary to the 
results obtained by Ryde & Svensson (2002). 

Next, we fit the pulses resulting from the assumed emis- 
sion in the expanding shell from i?min to iimax • In the Figure 
[7|we show the distribution of indices for all the sub-pulses 
that fitted to this pulse profile (again apart from the sub-jets 
seen at angles S^bs > 0.03). The distribution is now much 
broader, with certain enhancement at low values of n, but 
with no traces of n = 3 secondary peak, as found by Ryde 
& Svensson (2002). 



3.2.1 New sample from BA TSE 

Here we present the analysis of the new sample of BATSE 
data, in order to make further observational tests for the 
model. Now, contrary to Ryde & Svensson (2002), we study 
only the short duration bursts. The sample of bursts is listed 
in Tableland contains 83 GRBs, of durations from 0.067 
s to 1.954 s. The data from the energy band 50 - 300 keV 
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Figure 6. The distribution of the index n, fitted to the sub-pulses 
profiles due to the instantaneous emission at r = ro. 



were used and the bursts were selected so that the signal to 
noise ratio in this band was larger than 10. 

We fit the pulse profiles to the shape defined by the 
equation 1201 The data are in units of [cts/s] and we do 
not convert them to the photon flux, due to the large er- 
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Table 1. Sample of the short GRBs observed by BATSE. The bursts were selected on the basis of their signal to noise ratio, S/N > 10, 
in the energy band 50-300 keV. The total sample contained 182 bursts; the table lists only these bursts that fitted to the defined profile 
shape CEg.dni. 



ror bars in case of short bursts. The conversion to physical 
units would dramatically increase the uncertainty. 

In the Figure |H| we show the distribution of the index n 
in the observed sample. The plot contains only these bursts, 
that fitted to the defined profile shape. Only less than half 
of the originally selected bursts remained (83 out of 182) . 

As the Figure shows, no significant excess in the distri- 
bution is seen around the index n = 3. However, an excess at 
low values of n is seen quite clearly. The figure shows overall 
similarity to the theoretical prediction based on thick shell 
approximation while it is definitively different from Figure|S| 



3.3 Statistical distribution of bursts 

Below we show the results of the Monte Carlo calculations 
of the statistical distribution of GRB durations. In all the 
simulations we take A^'obs = 1000 observers (or, equivalently, 
1000 GRBs at different positions on the sky with respect 
to the observer). Each of the bursts consists of A'sub = 350 
sub-jets, which have the opening angles of AS^ub = 0.02 and 
Lorentz factor 7 — 100. The departure time of each sub-jet, 
tdep, is chosen randomly within the total activity time of 
the central engine, i.e. between and T^ut = 30 sec. The 
sub-jets have a quasi-Gaussian distribution within the jet of 
the total opening angle A^tot = 0.2. 

We compare the results for two cases: the instantaneous 
emission form a thin shell at r — Ro and the emission from 
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Figure 7. The distribution of the index n, fitted to the sub- 
pulses profiles due to the emission from the expanding shell from 
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Figure 8. The distribution of the index n, fitted to the profiles 
of the short GRBs observed by BATSE. 



the thick expanding shell, from r — Rmin to r = -Rmax- The 
infinitesimally thin emitting shell is located at ro = 10^^ cm, 
and the thick shell emits the radiation from i?min = 10^^ 
cm, to Rmax — 2 X 10^* cm. Below, in the Section. 13.3. II we 
discuss the distribution limited to these bursts that contain 
only the sub-jets, and in Section r3.3.2l we show the results for 
the bursts that include the emission from every sub-pulse. 



Figure 9. The simulation with N^iyg = 1000 random observers, 
and the jet consisting of A^sub = 350 sub-jets. The histogram 
shows the number of bursts (observer's positions) in the function 
of the total burst duration time. Each burst includes only the sub- 
jets caught on the line of sight. The emission is produced either 
by an extended shell from Rmin = 10^* cm, to -Rmax = 2 X 10^* 
cm (solid line) or by the thin shell at ro = 10^* cm (dashed line). 



seen either on axis or off-axis. The criterion for counting 
these sub-bursts as a part of the GRB is their flux above a 
background for the duration time Tgo, as observed by the 
BATSE detector. 



3.3.1 Statistics for nearly- on- axis observers 

As a first step, we will study the distribution of burst dura- 
tions under the (restrictive) assumption that only the sub- 
jets that are on the line of sight contribute to the GRB 
event. In other words, we neglect any emission that is ob- 
served off-axis, which could possibly contribute to gamma 
rays, assuming that those events would only give rise to 
the X-ray Flashes. We require that all the sub-jets that are 
contributing to the GRBs will satisfy the condition |H| with 
respect to the observer axis. 

The observer may catch a sub-jet on his line of sight 
either while looking very close to the main jet axis, with 
Sobs ^ (multiple sub-jets are seen), or while looking close 
to the edge. In the first case, the number of detected sub- 
jets is about 300, and in the latter case, the observer is 
most likely to catch one sub-jet (about 10% of observers). 
Obviously, the largest number of observers do not catch any 
single sub-jet. 

In the Figure]^ we show the histogram of burst dura- 
tions, for this simulation. The distribution has two peaks, 
located around 4 and 50 seconds, or 8 and 65 seconds, for 
an instantaneous or extended emission from the sub-jets, re- 
spectively. The exact location of the peaks depend on the 
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model parameters, and in particular the short GRB peak 
is sensitive to the value of the radius of emitting shell, as 
well as the opening angle of a sub-jet (cf. Eg. I12|l . For ex- 
ample, A^sub — 0.015 shifts the short burst duration peak 
to Tgo — 2.5 s, while the radius of ro — 10^'^ cm shifts this 
peak to Tgo = 0.4 s. The long burst peak does not depend 
on these parameters, but rather on the duration time of the 
central engine activity. If we take T^ut ~ 0.3 s, instead of 30 
s, the peak for long bursts is shifted to Tgo = 8 s. However, 
the peak for short GRBs does not change and is located at 
~ 4 s, so the shorter activity time results in the essentially 
single-peak distribution for only short bursts. 



3.3.2 Statistics for the off- axis emission included 

In this section we discuss the effect of the off-axis emis- 
sion that contributes to gamma rays. Since the observed flux 
drops very rapidly with the increasing inclination of the line 
of sight, as soon as the observer axis is located outside the 
sub-jet cone, it is crucial to determine the condition for the 
GRB to be observed. This is connected mainly with the flux 
limit of the detector. 

For example, BATSE will generate a burst trigger if 
the count rate in two or more detectors exceeds a threshold 
specified in units of standard deviations above background 
(nominally 5.5). The actual detection depends only slightly 
on the spectral hardness 

In the Figure [THl we show the results of the simulation. 
The off-axis emission of the sub-jets contributes to GRB 
whenever 



(21) 



where the background flux is arbitrarily chosen to be Fb = 
1.0 

The shortest GRBs are mainly due to the bursts that 
include only the off-axis emission. This is because the weak 
off-axis sub-bursts in general tend to have longer dura- 
tions, which sometimes prevents them from falling below 
the threshold. On the other hand, the weakest of these sub- 
bursts have too small flux at 100 keV, and are missed any- 
way. Therefore the total number of the sub-bursts that con- 
tribute to the left wing of this histogram is small and the 
only sub-bursts that are counted here are the off-axis sub- 
bursts. 

The sub-bursts caught on-axis (or nearly on-axis) give 
the main contribution to the longest GRBs. Since their emit- 
ted flux is always large, the threshold limit does not affect 
them. These sub-bursts are seen mostly at very small ob- 
serving angle Sobs, and a large number of them is detected. 
This means a large spread between the starting time of the 
first sub-jet and ending time of the last one, which translates 
into a long total event. 

The shape of the histogram is clearly not bimodal. A 
small excess of short bursts around T ~ 10 s appears only for 
the model with instantaneous emission at r = ro, however 
this feature disappears already in the model with expanding 
shell emission from rmin to rmax- The bimodal distribution is 
found by Toma et al. (2005) in the model with only instan- 
taneous emission at r = ro. However, these authors consider 
the case when the events are classified by the detector as 
a GRB on the basis of their spectral hardness. They select 
only these events for which the hardness ratio, i.e. the ratio 



200 




Figure 10. The results of the simulation with the emission from 
ofT-axis sub-jets included. The bursts include only these ofE-axis 
sub-bursts, for which F^^^y ^ 5-5x ■\/~Fg/Tgo. Other parameters 
are the same as in Fig. |2l 



of fluency defined as 5*1,1,^2 

= ""^ XT' F^dudT, is smaller 
than S'2kcV,30kcv/'S'3okcV,4ookcV ~ 0.3. In this case the con- 
tribution from the short bursts is much more pronounced. 



3.3.3 Distribution of redshift 

The distribution of the long bursts with redshift is now 
thought to trace the star formation in an unbiased way 
(Jakobsson et al. 2005). Therefore the probability of find- 
ing a bursts with a redshift z is given by: 



dP{z) = AD^R(z)dz 



dr 
dz 



(22) 



where R[z) is the star formation rate in the comoving frame 
(Porciani & Madau 2001), D is the cosmological distance 
(Chdorowski 2005), A is the normalization constant and 
dr/dz = c/Ho/ ^VLmi'i + z)'* + We adopt the cosmo- 
logical parameters of Om = 0.3, f^A = 0.7 and Ha — 65. The 
normalization is calculated from the condition that: 



P{z)dz - 



(23) 



Here we check the infiuence of the redshift distribution 
on our model. For each of the 1000 bursts we randomly 
choose its redshift, within a distribution with a cut-off at 
Zmax = 4. The Figure 1111 shows the resulting statistics of 
bursts. 

The duration distribution is now more smeared than 
in case of a constant redshift, z — 1.0, and shifted toward 
longer durations. The total number of bursts that satisfy 
the condition 1211 1 is similar to that in Fig. 1101 and about 
~ 700 out of a 1000 bursts exhibit a large fiux in gamma 
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Figure 11. The same as in Fig. IIUI but with a random redshift 
distribution included. The redshift for each burst was chosen be- 
tween and 4, with a probability distribution that traces the star 
formation rate, Rsvi^yi = 0.3, Ha = 0.7). 



rays. However, the durations of the bursts with the distri- 
bution of redshifts tend to be more equal, due to the linear 
dependence on 1 -I- z. This makes the peak of the histogram 
flatter. The shift toward longer durations results from the 
fact the redshift distribution traces the star formation rate 
according to the function -Rs/2 in Porciani & Madau (2001) 
so the number of events with z < 1.0 is very small while sat- 
urates for large z. Therefore the long bursts can be longer, 
due to the contribution of bursts with z > 1.0, while the 
short bursts have similar durations to these in Fig. IIUI 



4 DISCUSSION AND CONCLUSIONS 

The scenario in which the central engine of a GRB operates 
like a "shot-gun" emitting multiple narrow bullets was origi- 
nally suggested by Heinz & Begelman (1999). Also, Fenimore 
et al. (1999) showed that the variability of GRB may be due 
to the breaking of the local spherical symmetry of the emit- 
ting shells. In this case the duration of the total event was 
determined by the time of activity of the central engine, and 
therefore connected with the progenitor type. The observa- 
tional evidence for the association of the long GRBs with 
supernovae (Bloom et al. 2002; Stanek et al. 2003) support 
the coUapsar model. The simulations of the jet propagat- 
ing through the envelope of the collapsing star (e.g. Aloy 
et al. 2000; MacFadyen et al. 2001) show that the jet has 
an ultimate Lorentz factor of about 100-150 and an open- 
ing angle wider than ~ 10° . The jet does not have a uniform 
structure in its density, velocity, pressure etc. These physical 
quantities have a certain distribution along with the angu- 
lar distance from the jet axis, however, the gradients are not 
very sharp, and it is not certain if this scenario could result 



in a separate "bullets" released from the center during the 
time of activity. 

One of the possible scenarios is the quasi-universal 
structured jet model (e.g. Rossi et al. 2002; Zhang & 
Meszaros 2002; Lloyd-Ronning et al. 2004), in which the ob- 
served diversity of GRB properties is explained by a different 
viewing angle. The jet has a standard geometrical configu- 
ration, with a stable outflow channel but jet properties (e.g. 
Lorentz factor, energy flux) are distributed non-uniformly 
per solid angle within the cone. 

This is in contrast with the simple jet models, in 
which the energy distribution within the cone was uniform 
(Rhoads, 1997; Panaitescu et al. 1998). These simple jets 
have been used to successfully explain the breaks in the GRB 
optical afterglow lightcurves (Moderski et al. 2000; Huang et 
al. 2000; Granot et al. 2002). Since the energy is coUimated 
to different degrees in various bursts, the jet opening angle 
is the main parameter that drives the jet properties. Which 
mechanism drives the various coUimation degree in various 
GRBs remains to be explained. 

The model studied in the present work invokes a kind 
of a "shot-gun" model, with a fixed opening angle, and a 
Gaussian distribution of the subjet probability around the 
symmetry axis. In the case of assumed instantaneous emis- 
sion the model is effectively quite close to a "shot-gun" sce- 
nario while in the case of extended emission (i.e. thick shell 
approximation) the emitting gas has definitively a jet-like 
appearance. However, the model is different from the sin- 
gle structured jet since the distribution of the directions of 
sub-jets contains the element of randomness. The underly- 
ing assumption is that the process of the jet formation is 
not stable, and an open jet channel may close up and be 
replaced with a new one. In this model the GRB properties 
depend significantly on the observer's position. 

We first analyzed the profiles of single pulses coming 
from separate sub-jets. 

Observationally, such profiles were studied by Ryde & 
Svensson (2002) for Tgo > 2 s bursts in BATSE sample. 
They found a good representation of single subpulse profiles 
with equation 120L and noticed a two-peak distribution of 
the index n: a higher peak around n = 1 and a lower peak 
around n = 3. We completed their analysis for T90 < 2 and 
found one peak below n = 1 but no significant traces of the 
second one. If we restricted ourselves to the data with the 
highest S/N ratio (above 30) perhaps we could see a slight 
relative enhancement at n ~ 3, but we have only 11 such 
events in the sample. 

We compared these results with the theoretical profiles. 
Assuming instantaneous emission we obtained values of n 
between n = 4 andn = 6 . However, the emission from 
thick shell with 7?min = lO^** cm and iimax — 2 x 10^* cm 
results in pulse profiles that are slightly wider in their rising 
part and have a different decay profiles. The corresponding 
distribution of the index n is much broader and peaks below 
n = 1. Therefore, the thick shell approximation seems to be 
a better representation of the observational data within the 
frame of the adopted model. 

We next analyzed the question whether the adopted 
model can explain the apparent existence of the two classes 
of bursts: short and long. 

We found that the multiple sub-jet model can repro- 
duce the bimodal distribution of GRB durations only for 
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some specific conditions and parameters. In particular, this 
kind of distribution is reproduced if only the sub-jets that 
are seen on the line of sight can contribute to the GRB. 
When this assumption is released and all the sub-jets can 
contribute to the GRB (provided they give a sufficient flux in 
gamma-rays), the bimodal distribution is smeared. Further- 
more, the bimodal distribution vanishes completely when- 
ever we assume the emission from a thick shell from -Rmin 
to -Rmax instead of instantaneous emission at _Ro 

In the sub-jet model studied in this paper the duration 
time of the GRB depends on two parameters: (i) the radius 
of the emitting photosphere of a shell, ro (or, more specifi- 
cally, its combination with the sub-jet opening angle, ASgub) 
and (ii) on the time of activity of the central engine, Tdur. 
The case (i) refers to the single sub-pulses detected by the 
observer, while the case (ii) is for the multiple sub pulses. 

An important ingredient of the model is the multiplicity 
of the sub-jets. This is parameterized by the ratio of the solid 
angle covered by all the sub-jets to the solid angle covered 
by the whole jet: ^ = (iVsub x M^^^) / M^^^. If this ratio 
is larger than unity, i.e. we have multiple sub-jets scattered 
with some, say, Gaussian, distribution around the jet axis, 
most of the bursts contain multiple pulses and the chance of 
detecting a single pulse is large only very close to the edge of 
the jet. For the standard parameters used throughout this 
article we had ^ = 350 x (0.02)^/(0.2)^ = 3.5 The duration of 
the GRB is governed then by the activity time Tdm, and the 
simulations show that the bimodal distribution is possible 
only for two distinct values of this activity time. 

We have a different situation for ^ <^ 1. Most of the 
bursts will then contain only single pulses, and their dura- 
tion is governed by the photospheric radius ro. We checked, 
that for extremely low multiplicity (A^sub = 10, which gives 
^ = 0.1) the bimodal distribution duration can be recovered, 
both for instantaneous and extended shell emission (see Fig- 
ure However, a question arises, if in such a case we can 
really talk about the "multiple sub-jet model". Also, the 
observed GRB pulse profiles would look totally different, as 
we would have mostly single or double spikes even for ob- 
servers looking very close to the jet axis. The observations of 
long GRBs show something different: many of them exhibit 
a narrow substructure rather than couple of spikes. 

To sum up, we suggest that an explanation of the bi- 
modal duration distribution of the observed GRBs is possi- 
ble rather on the basis of two independent kinds of sources 
responsible for short and long events. This is further justi- 
fied by the fact that apart from many similarities, there are 
also systematic differences between long and short bursts, 
such as their spectral hardness or the systematically lower 
fluence in the case of short bursts. 
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Figure 12. Histograms for models with very small sub-jet multi- 
plicity: A^sub = 10, A^sub = 0.02, A^tot = 0.2. A: emission from 
ro = 10^'' cm (dashed line) and from -Rmin = lO'^* cm to -Rmax = 
2 X 10^'* cm (solid line); B: emission from -Rmin = 1-66 X lO'^^ cm 
to -Rmax = 3.33 X 10^3 cm. 
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